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Abstract 



OO 



We extract the K + — ► 7r + e + e~ amplitude scale at q 2 = from the recent Brookhaven E865 
high-statistics data. We find that the q 2 = scale is fitted in excellent agreement with the 
theoretical long-distance amplitude. Lastly, we find that the observed q 2 shape is explained 
by the combined effect of the pion and kaon form-factor vector-meson-dominance p, lo and 



ly-) ' 4> poles, and a charged pion loop coupled to a virtual photons e + e~ transition, 

m 1 PACS numbers: 13. 20.Eb and 13.25.Es 

o 
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^ ■ I. INTRODUCTION 

a: 

The decay K + — > 7i + e + e has been studied theoretically for many years. There is general agree- 



ment that the process is dominated by the long-distance "bremsstrahlung" graphs, in which a 
virtual photon is radiated by the charged kaon or the pion, away from the strangeness-changing 
weak vertex. However, qualitative agreement with the available experimental evidence on the 
branching ratio and the form factor (i.e. the q 2 dependence) has not previously been achieved. 
Here, the virtual photon mass, q 2 , is given by (p + + p~) 2 , where p± are the four-momenta of 
the e . Also, most calculations have several unknown parameters, making firm predictions dif- 
ficult. Experimentally, the situation has been dramatically improved by the recent availability, 
from Brookhaven experiment E865 IJ, of accurate data on both the branching ratio and the q 2 
dependence. In this paper, we fit these data to about 5% in a calculation based on a specific 
model with no free parameters (except for relative signs), which fits other related data. 

II. AMPLITUDE AT q 2 = 

The E865 experiment 0] at Brookhaven had measured a new value of the branching ratio for K + — > 

7r + e + e~. The value depends on the form of the extrapolation to q 2 = 0. We choose to extrapolate 

with a form factor quadratic in q 2 , which gives BR(K + — > ir + e + e~) = (2.99 ± 0.06) x 10 -7 (see 

Table 1, column 2 of Ref. [|TJ). From this, one extracts the invariant amplitude defined by 
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giving 



|A exp (0)| = (4.00 ±0.18) x 10- 9 Ge\r 2 
found from the 3-body phase-space integral for the rate 



(2) 



T(K + -> vr+e+e" 



0.2020 m x\ A ^}\ 2 = (1.59 ± 0.03) x 10- 23 GeV . 
3(47r)' :i 



(3) 



As shown in |T[], a quadratic fit provides a better agreement than a linear one for the same data. 
We will return to the q 2 dependence later. 
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Figure 1: Graphs for K + — > 7r + e + e~ through a virtual photon, (a) and (b) are long-distance 
graphs, (c) is a short- distance graph and (d) is a pion loop term. In each graph, the blob denotes 
the weak (strangeness-changing) vertex and the wavy line is an off-shell photon. 



The long-distance (LD) chiral low-energy model, used to calculate the bremsstrahlung graphs 
of Fig. 1(a) and (b), predicts [0 



A 
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(7T+ | H w | K+) 



m K+ ~ 



F^+{q 2 ) - F K+ (q 2 
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where F 1T +{q 2 ) and F K +(q 2 ) are the pion and kaon electromagnetic form factors. At q 2 = 0, Eq. 
becomes 
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We now present the evaluation of the matrix element in Eq. ([5]), describing how we obtain 3.9 x 
10~ 9 GeV~ 2 for Ale>(0), reasonably close to the experimental result in Eq. (0). The matrix element 
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I (tt + I Hy/ | K + ) | is well established; in Ref. [^J its value was deduced theoretically and was 
confirmed by comparison with values from ten measured kaon decays, Kg — > 2tt°, K — ► 37r, 
Kl s — ► 27 and Kl — > 7r°27, all of which are consistent with | (n + \ H\y \ K + ) \ = \ (tt° \ Hy/ \ 
K L ) |^ 3.5 x l(r 8 GeV 2 . Specifically the K s -> 2tt° rate T gives || 

I (2tt° I I K s ) |= m xv /l6vrr/p = (37.1 ± 0.2) x 10" 8 GeV , (6) 

where p is the three-momentum of a ir° in the Ks rest frame. Using partially conserved axial 
currents (PCAC), this predicts for the pion decay constant f„ ~ 93 MeV the LD scale 

I (tt+ I H w I K+) I w I / t (2tt I if w I K s ) I « 3.5 x 10- 8 GeV 2 . (7) 

For F n +(q 2 ), there are many experimental measurements for both positive and negative g 2 , 
including our region < q 2 < 0.125GeV 2 . The experimental data close to the region of relevance 
here are shown in Fig. 2. As expected from Vector Meson Dominance (VMD), these data are well 
described by a p pole; the curve in Fig. 2 is the p-pole, i.e. VMD prediction, 

The pion charge radius r n = ^ < r 2 >, where r 2 = 6dF(q 2 )/dq 2 \ q 2 =0 , which in VMD equals 
0.628 fm, is in agreement with the experimental value || r T = (0.63 ±0.01) fm. Since both F n (q 2 ) 
and 7V agree well with the data, the p-pole expression can either be regarded as an input to our 
calculation or as an interpolating function between the experimental points. Note that up to a 
small G-parity violation, u and (ft do not contribute to F n +. 
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Figure 2: Pion electromagnetic form-factor data ^ and the prediction of the p pole from VMD, 
Eq. ([§[). The p width is included here; it has practically no effect in the range of q 2 relevant for 
K + — > 7r + e + e~. 
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For the kaon case, no such comparable data exist, so we need a model for Fx+{q 2 )- In view of 
the success of the VMD picture for F„-+(g 2 ), we use the same model for F K +(q 2 ), resulting in the 
p, u and (f> pole structure 



F K+ {q 2 




(9) 



where g pee = 5.03, g^ ee = 17.06, g^ ee = 13.24 (derived from the e + e~ decay widths), and the 
p°K + K~, uK + K~, (pK + K~ SU(3) coefficients are |, ~, -j=, respectively. We deduce the pKK 
coupling constant from g piT7T (obtained using the p width) via SU(3). We link the normalization 
F K +{0) = 1 to iV = 0.037GeV 2 in Eq. © leading to Grip^ = 0.94,3.08, 1.99. We then find that 
the charged K radius is r K + = 0.557 fm, in very good agreement with the data pj. The above 
lead to the prediction r 2 — r 2 K = 6 (dF^+fdq 2 — dF K + /dq 2 )\ q2=0 = 0.063 fm 2 , consistent with the 
experimental value |8| (0.100 ± 0.045) fm 2 . Substituting into Eq. (Q) the above values for the 
derivatives of F(q 2 ) at q 2 = (see Eqs. || and |9|), together with Eq. (0) predict 

A^ D (0)| = 3.9 x 10- 9 GeV- 2 . (10) 

Folding in the p width into the form factors and their derivatives has a negligible effect for the 
range of q 2 relevant to us here, i.e. q 2 between and 0.125GeV 2 . 

Our assumptions for F n +(q 2 ) and Fx+(q 2 ) are consistent with a specific model, which has the 
VMD and structure built in, namely the quark-level linear a model (LctM) as discussed in 
Refs. H, [l(J. The pion electromagnetic form factor, obtained from a u, d quark triangle graph 
for a 7r + probed by a photon (with pseudoscalar -wqq coupling (775), predicts 

nonperturbatively normalized to ^+(0) = 1 via the quark-loop version of the pion decay constant 
fn, combined with the quark-level Goldberger-Treiman relation f n g = m q . To extract the q 2 
dependence from Eq. (|11|) while integrating out the quark momentum p, one can differentiate 
Eq. QII]), using r 2 (g 2 ) = QdF n+ (q 2 ) / dq 2 for g 2 = (2rr) 2 /N c and r n+ = l/m q w 0.63 fm @, 0: 

rl+(q 2 ) = Q^dx x(l - x) [m 2 - x(l - x)q 2 }' 1 . (12) 

Performing the above integration analytically and then making a Taylor series expansion keeping 
only the leading terms with y = q 2 /m 2 , one finds 

+ | + + ^ + • (13) 
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Since m q pa 171^/3, so m 2 /m 2 ~ 6, one can approximately express Eq. ([13]) in VMD language for 
the form factor itself in the low q 2 region as in Eq. @. Lastly, the coupling g pn7T = 6.04 is taken 



from the measured p width, with the LaM correction to VMD predicting [0 g^/Opee — 6/5, 
very close to the data ratio 6.04/5.03. Then from SU(3) one gets g p x+K- — 3.02. A test of this 
approach is to extract g^x+K- — 4.58 from the measured partial width of <fi — > K + K~ Q, from 
which we find, using SU(3) that g P K+K- = 3.24. This value agrees to better than 10% with its 
value 3.02 above. 

In the work of Ref. 0, the short-distance (SD) term was estimated to cause about a 20% 



reduction in the LD term. But more recent work by Dib, Dunietz and Gilman |TTJ, based on a 
much heavier top-quark mass and also including QCD corrections, suggests that the SD term is 
much smaller. We therefore neglect the SD term of fig. 1(c) in this work. 

Our predicted scale, then, is the LD amplitude of Eq. (pi]), which agrees with the data, Eq. (|j), 
to about 5%. 

III. FORM FACTOR 

To extract the form factor for K + — > 7r + e + e _ from Eq. (f|), we need the q 2 dependence of the 
electromagnetic form factors, F n +(q 2 ) and F K +(q 2 ). In applying Eq. (Q), we take F n +(q 2 ) from the 
data in Fig. 2, parameterised using the p pole (Eq. (§)) and for F K +(q 2 ), we use the equivalent 
expression, Eq. @. We emphasise that we do not rely heavily on VMD for F n +(q 2 ); we use 
essentially the empirical values. The main reason to show the VMD fit to F 7r +(g 2 ) is to establish 
that this theory fits the data well, so it is reasonable to use it to obtain Fx+(q 2 ), to apply in 
Eq. (|). 

With F 7r +(g 2 ), F K +(q 2 ) and (ti + \ H w \ K + ) taken as above, the LD part, Eq. (§), is completely 
determined. The resulting form factor, Ald, rises with q 2 , in agreement with the data, but gives 
only about 30% of the observed rise in the form factor, | F | 2 , between q 2 = 0.03 and 0.12 GeV 2 . 
This is shown by the dashed line in Fig. 3. 



Another possible contribution has been discussed by Ecker et al. [|ll| and by D'Ambrosio et 



al. [0], namely the charged pion loop term, Fig. 1(d). This term is derived in Refs. [12| and 



H|] using dimensional regularisation. We use here the expression for W 77 ^ from Ref. |T3[ without 



any additional polynomial. Also, we do not include terms from Refs. Jl2j and [13] other than 



the loop term; these contributions are calculated explicitly in our amplitude A^. Evaluation of 
the pion loop term requires a knowledge of the K + — > 7r + 7r + 7r~ amplitude. In Refs. [0 and 



[ 13|1 this is taken from experiment. Ref. |14| shows that a current algebra-PCAC approach for 



K + 7r + 7r + 7r agrees with K 3n data within 5%, so the g 2 -dependent part of the loop in Refs. 



[ |12j l and [|13j is, in fact, compatible with the methods used here. The relative sign of the pion 
loop term and the dominant part of the amplitude, which may be considered as a parameter, was 
already established by the experiment PJ. 



5 



4.5 



4 
2 



3.5 



2.5 



1 .5 



0.5 




0.02 



0.04 



0.06 



0.08 



0.1 



0.12 



q 2 (GeV 2 



Figure 3: Form factor squared, \ F \ 2 , as a function of q 2 for K + — > 7i + e + e~ . The dashed curve 
shows A^d and the solid line shows A^d + A n \ oop . The black dots are the experimental data 
The theoretical curves are normalized to 1 at q 2 = . Complex masses are used here; inclusion of 
the imaginary parts of the masses has a small, yet visible effect on the high q 2 region in this figure 



Adding the amplitude, A^op, from the pion loop term to Ald gives the form factor shown 
as the solid line in Fig. 3. This agrees quite well with the data. The pion loop term makes a 
negligible contribution to the amplitude at q 2 = 0, so it does not disturb the agreement with 
experiment of the scale, A(0), discussed in Sec. II and resulting in Eq. (|TOl) . 

IV. CONCLUSIONS 

In summary, we have shown that the recent Brookhaven data for K + — > 7r + e + e~ are in qualitative 
agreement, both in the q 2 = scale and the form factor, with a calculation in which two processes 
are included. The experimental amplitude, |A exp (0)| = (4.00 ± 0.18) x 10~ 9 GeV~ 2 (Eq. ©), is 
about 5% around the simple VMD bremsstrahlung prediction of |Ald(0)| = 3.9 x 10~ 9 GeV~ 2 in 
Eq. (|T0|), having dropped the SD term according to Ref. fliT] . 

This agreement is as good as can be expected from such a simple model. In particular folding 



in a K* pole through K A 



7T + K 



*0 



7T 1 e e 



and a small, about 10% SD contribution to 



the amplitude []T5 , may modify our prediction slightly. Note that as more terms are added, one 
increases the danger of double counting. Each of the graphs discussed here has been calculated 
before, though all have never been taken together in one calculation. This level of agreement is 
reasonable in view of the simplicity and combination of the models used. 

There are no free parameters in our treatment, except for the relative sign between the pion 
loop and the rest of the amplitude and a relative sign between the couplings in Eq. @; all 
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quantities are taken from experiment, or from models that are known to fit other low-energy kaon 
or pion processes well. We do not attempt to estimate theoretical errors on our calculation, since 
these are to some extent dependent on one's point of view. For example, the pion electromagnetic 
form factor, shown in Fig. 2, could equally well be regarded as taken from experiment, with 
its associated errors, or from the VMD prediction which is essentially free from error within the 
model. Similarly, the value of | (tt + | Hy/ \ K + ) | is derived from theoretical arguments, but is 
confirmed by the fit to no less than ten other kaon decays [2,3]. Either of these sources could be 
regarded as the origin of our value for | (tt + | Hw \ K + ) | and hence as the source of error in it. 

Concerning the recent paper [|16j on K + — > 7i + ee, we agree with its recommended form factor 
structure (as in our Eq. ([|); compare with its Eqs. (2), (7) and (8)). However we disagree about 



the scale of Eqs. (7) and (8) in |16| which seems more akin to the small AI = 3/2 K — > -kit 
tree level transition rather than to the much larger K° — ■> tttt AI = 1/2 LD transition, as in our 
Eqs. (g) and (0). 

As for Chiral Perturbation Theory (ChPT) applied to K + — > 7r + ee |17[| , the required 0(p 6 ) 
expression contains about 100 parameters. Moreover, the underlying A J = 1/2 LD scale as in our 
Eq. (0) disappears for Nc — > oo |17| . 

ACKNOWLEDGEMENTS 

We are grateful to T. Goldman, G.J. Stephenson Jr., P. Singer and F. Krauss for useful 
discussions. We acknowledge support from the US DOE. GE would like to thank the US-Israeli 
Binational Science Foundation and the Fund for the Promotion of Research at the Technion for 
partial support. MDS appreciates partial support from the Technion. 



References 

[1] E865 collaboration, R. Appel et ai, Phys. Rev. Lett. 83, 4482 (1999). 
[2] G. Eilam and M.D. Scadron, Phys. Rev. D31, 2263 (1985). 

[3] M.D. Scadron and R.E. Karlsen, Intersections 91 Conference, AIP Conference Proceedings 
no. 243, ed. W.T.H. Van Oers, 1992, p. 627; Nuovo Cimento A106, 113 (1993). 

[4] Particle Data Group, D.E. Groom et al, Eur. Phys. J. C15, 1 (2000). 

[5] R.E. Karlsen and M.D. Scadron, Phys. Rev. D44, 2192 (1991). 

[6] S.F. Bereshnev et al, Sov. J. Nucl. Phys. 18, 53 (1974); 24, 591 (1976); A. Quezner et al, 
Phys. Lett. B76, 512 (1978); LB. Vasserman et al, Sov. J. Nucl. Phys. 33, 368 (1981); E.B. 



7 



Dally et al, Phys. Rev. Lett. 48, 375 (1982); A.F. Grashin and W.V. Lepeshkin, Phys. Lett. 
B146, 11 (1984); L.M. Barkov et al, Nucl. Phys. B256, 365 (1985); S.R. Amendolia et al, 
Nucl. Phys. B277, 168 (1986). The last of these references contains an excellent summary of 
all earlier work. 

[7] T.H. Bauer, R.D. Spital, D.R. Yennie and F.M. Pipkin, Rev. Mod. Phys. 50, 261 (1978); 
Erratum-z&zd. 51, 407 (1979); L. Ametller, C. Ayala and A. Bramon, Phys. Rev. D24, 233 
(1981); C. Ayala and A. Bramon, Europhys. Lett. 4, 777 (1987). Eq. (11) in the latter 
reference has VMD forms similar to our results in Eqs. (|| |9|) above. 

[8] S.R. Amendolia et al, Phys. Lett. B178, 435 (1986). 

[9] R. Delbourgo and M.D. Scadron, Mod. Phys. Lett. A10, 251(1995); T. Hakioglu and M.D. 
Scadron, Phy. Rev. D43, 2439 (1991). 



[10] A. Bramon, Riazuddin and M.D. Scadron, J. Phys. G 24, 1 (1998); |hep-ph/97092T4] also 



see N. Paver and M.D. Scadron, Nuovo Cimento A78, 159 (1983); The equivalent result 
= N c /47r 2 f% was earlier obtained by R. Tarrach, Z. Phys. C2, 221 (1979); S.B. Gerasimov, 
Sov. J. Nucl. Phys. 29, 259 (1979); V. Bernard, B. Hillar and W. Weise, Phys. Lett. B205, 
16 (1988). 

[11] CO. Dib, I. Dunietz and F.J. Gilman, Phys. Rev. D39, 2639 (1989), based on the formalism 
of F.J. Gilman and M. Wise, Phys. Rev. D21, 3150 (1980); See also L. Bergstrom and P. 
Singer, Phys. Rev. D43, 1568 (1991). 

[12] G. Ecker, A. Pich and E. deRafael, Nucl. Phys. B291, 692 (1987). 

[13] G. D'Ambrosio, G. Ecker, G. Isidori and J. Portoles, JHEP, 8, 4 (1998); G. Ecker, private 
communication. 

[14] R.E. Karlsen and M.D. Scadron, Phys. Rev. D45, 4108 (1992). 

[15] See L. Bergstrom and P. Singer in [IT] and P. Singer, |hep-ph/960742!i . 



[16] P. Lichard, Phys. Rev. D61, 073010 (2000). 
[17] For a summary see: G. Isidori, |hep-ph/001iOT7| . 



8 



